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The transmission of digital data over telephone channels has been con- 
sidered previously in the literature, and the effects of Gaussian noise have 
been analyzed. With experience, however, it has become apparent that while 
a background of Gaussian noise is present, the limiting noise is not Gaussian 
but impulsive in nature. It consists of bursts of high amplitude which occur 
at random considerably more often than is predicted by the rms value of the 
Gaussian background noise. 

Measurements of the statistics of this noise have been initiated, and some 
results have been reported. In this paper, a model for the noise is constructed 
to be reasonably consistent with these measurements without becoming too 
complex to be handled analytically. Various modulation systems are ana- 
lyzed to determine their performance in such a noise environment. Condi- 
tional error rates, in terms of the average number of bit errors per noise 
burst, are determined as functions of a convenient signal-to-noise ratio 
which is defined. The systems are ranked as to their performance in such a 
noise environment, and the ranking is found to be the same as that for 
Gaussian noise. 

The improvement to be gained by employing complementary delay net- 
works is investigated. Networks with linear and sinusoidal delay charac- 
teristics are considered. 

I. INTRODUCTION 

The limiting noise with regard to transmitting digital data is impulsive 
in nature. Consisting of bursts of high amplitude, it occurs at random 
considerably more often than is predicted by the rms value of the 
Gaussian background noise. In this paper, an analysis is made of digital 
transmission, by various modulation procedures, in such a noise environ- 
ment. 



* Taken from the dissertation submitted to the Faculty of the Polytechnic 
Institute of Brooklyn in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy, 1964. 
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The block diagram of a generalized data transmission system includ- 
ing a telephone channel is shown in Fig. 1. The data signal consists of a 
train of either ideal impulses or square pulses, each pulse having either 
positive or negative polarity depending on whether it represents a mark 
or a space. The transmitter consists of a low-pass filter, which band 
limits the data signal, followed by a modulator and a band-pass filter. 
The low-pass filter is required to prevent "foldover" distortion in modu- 
lation. The band-pass filter restricts the transmitted signal to the range 
of frequencies passed by the channel. It avoids the waste of transmitted 
power in signal components which will not be received and also includes 
the channel splitting filters which prevent crosstalk into frequencies re- 
served for other channels. The low-pass and band-pass filters also shape 
the signal, giving it the form desired for transmission. The transmitted 
signal is applied to the channel, where it is contaminated by additive 
noise. The combined signal and noise enters a receiver which consists 
of a band-pass filter and demodulator which is generally followed by a 
low-pass filter and a synchronous decision device. The band-pass filter 
removes components of the noise outside the band of the signal, in addi- 
tion to shaping the received signal. The low-pass filter, which is not 
required in every system, removes the demodulation products which lie 
outside the band of interest, and may also provide further signal shap- 
ing. The decision device samples the combined signal and noise at its 
input at discrete sampling instants, and on the basis of each sample, 
produces a mark or space symbol. 

The channel is band limited, and may have an amplitude versus fre- 
quency characteristic which is not fiat and a phase versus frequency 
characteristic which is not linear. However, for the purposes of analysis, 
the exact characteristics of the channel may be lumped with those of 
the receiver band-pass filter. The channel may then be considered as 
having the characteristics of an ideal band-pass filter, with unity gain 
and zero phase shift in the band of interest. 
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Fig. 1 — Generalized data transmission system. 
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The very broad objective, in considering such a problem, is to design 
the transmitter and receiver so that the number of errors caused by the 
noise is minimized, under the constraint that the average power of the 
transmitted signal is limited. A corollary objective, and obvious pre- 
requisite, is to derive an expression for this error rate as a function of 
some readily measured signal-to-noise ratio, with the system charac- 
teristics as parameters. Then, one or two simple measurements of the 
noise on the channel will suffice to predict the error rate which can be 
expected of any system. 

In this paper, various existing modulation systems are considered. 
The specific problem, then, is to find a suitable mathematical model for 
the noise, and for each system, to determine the response of the demodu- 
lation procedure to the combination of the signal and such noise. The 
processed noise at the output of the demodulator and filters is the source 
of the errors by the decision device. The optimization procedure con- 
sists of finding the decision criterion and filter characteristics at the 
receiver which minimize the number of errors caused by the noise, under 
the constraint of average power limiting of the transmitted signal. For 
each modulation system, this procedure specifies the receiver design, 
except for the phase characteristic of the receiving filter. In order to 
avoid intersymbol interference at the decision device, a specific func- 
tional form for the signal at that point is required. Given the filter char- 
acteristics at the receiver, as determined by the optimization procedure, 
the filter characteristics at the transmitter are made such that the re- 
sponse of the over-all transmission path to the applied data signal is the 
specific signal required at the decision device. This specifies the trans- 
mitter design for each modulation system. 

In this paper, the following modulation systems are analyzed to de- 
termine their performance in the presence of impulsive noise: 

(1.) Double sideband AM with coherent detection 

(2.) Single sideband AM with coherent detection 

(3.) AM with envelope detection 

(4.) Frequency shift keying with frequency discrimination 

(5.) Binary phase shift keying with differentially coherent detection 

(6.) Quaternary phase shift keying with differentially coherent de- 
tection. 
Conditional error rates, in terms of the average number of bit errors per 
noise burst, are determined as functions of a signal-to-noise ratio which 
is defined in the following section. The systems are ranked as to their 
performance, in order of increasing conditional error rate. This ranking 
is found to be the same as that for performance in the presence of 
Gaussian noise, and agrees with the ranking experienced in actual usage. 



1702 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1905 

In the analyses, it is found that the phase characteristics of the re- 
ceiving filter affect the performance of the system. This phenomenon! , 
which is not present when the noise is Gaussian, leads to the considera- 
tion of complementary delay filters. Three types of delay filters are con- 
sidered, and the improvement which results from their use is found, as a 
function of the maximum delay of the filter, for each of the modulation 
systems. 

II. MODEL FOR THE IMPULSE NOISE 

As previously described the various data transmission systems include 
modulation and demodulation of a carrier frequency. The data receivers 
each include a band-pass filter. For the purposes of analysis, the actual 
characteristics of the channel are lumped with those of the band-pass 
filter, and the channel is considered as having the characteristics of an 
ideal band-pass filter. With double-sideband transmission, these charac- 
teristics are made symmetrical about the carrier frequency. The impulse 
noise present on such a channel consists of bursts of carrier frequency 
w c , with random phase ^, occurring randomly in time. When single- 
sideband transmission is used, the characteristics are asymmetrical and 
the Hilbert transform of the envelope is also present. This however, 
merely modifies the envelope and phase of the noise burst. The actual 
noise, as it occurs in the telephone plant, is wideband, with a spectrum 
covering many channels. The noise burst at the output of any one chan- 
nel is the response of the channel to the wideband noise and contains 
only a small portion of the total spectrum. It is reasonable to assume 
therefore, that the spectrum of the noise burst at the output of the 
channel is essentially determined by the channel characteristics, with 
the original spectrum of the wideband noise having little effect. Under 
this assumption the envelope of each noise burst is the same, except for 
a random amplitude K. This is of course, an approximation. The spectra 
of the noise bursts do vary somewhat. However, in order to be useful, a 
model must be reasonably consistent with the observed phenomena 
without becoming too complex to be handled analytically. Approximat- 
ing the envelopes of the noise burst by a representative time function 
n (t) is a very reasonable approximation. At this point in the analysis, 
the envelope n (t) is not limited to any particular function. A general 
expression for the error rate is derived, into which any specific function 
may be inserted. Then in order to obtain numerical results, a specific 
function is assumed. Up to that point, however, the analysis is quite 
general. 
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A single noise burst of carrier frequency n(l), occurring at time t = t, 
may be represented in the form: 

n(t) = Kn (t - t) cos ( Ue t + $) (1) 

where K, t, and \j/ are three independent random variables associated 
with each burst. For ease of notation, the envelope n (/) is normalized 
to have an energy of two watt-seconds. Then, the energy e of the noise 
burst n(t) is equal to K 2 . 

Measurements have been made on representative telephone channels, 
to determine the statistical characteristics of the noise. 1 - 3 As a result, 
there have been functional forms for: 

F(x) = Pr[K > x] (2) 

suggested by Fennick, 3 Mertz, ' and others, which fit the measured data 
reasonably well for the range of K large enough to cause errors. These 
are families of functions F a (x), where the parameter a may vary from 
channel to channel but is constant for any one channel. One such family 
of functions, suggested by Mertz, which fits the data reported on by 
Fennick, is computationally suited to this development. A reference 
amplitude Ko is chosen sufficiently low in value that any noise burst 
which might cause an error would have an amplitude greater than K . 
The relative frequencies of amplitudes greater than K are measured and 
plotted. That is to say, the conditional probability 

?r[K > x | K > Ko] 

is plotted versus the relative amplitude 

20 \og 10 (x/K )db. 

These plots are straight lines on logarithmic paper, as shown in Fig. 2. 
The negative slope of each line, in decades per 10 db, is the parameter 
a. These plots satisfy the equation 

Pr[K > x | K > Ko] = (K /x) 2a . (3) 

When describing the impulse noise, it is sufficient to consider only those 
bursts with amplitude greater than K ; if Ko is chosen low enough that 
bursts with smaller amplitude do not cause errors, these smaller bursts 
may be ignored and artificially assumed not to occur. Then, the family 
of functions F«(a) is given by: 

F a (x) = Pr[/v > x] = (K /x) 2a for x > K . (4) 

The majority of the measured distributions have values of a between 1 
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Fig. 2 — Distribution of noise burst amplitudes. 

and 2.5, with occasional higher values. High values of a correspond to 
channels with a small percentage of high amplitude bursts, and thus 
result in low conditional error rates. Lower values of a correspond to 
channels with a greater percentage of high amplitude bursts, resulting 
in higher conditional error rates. A distribution with a equal to 1 would 
have an infinite variance; the mean of K 2 , which is the average energy 
in a noise burst, would be infinite. A value of 1 is a lower bound on a, 
resulting in an upper bound on the conditional error rate. 

A function which is more useful in the ensuing analysis is the distribu- 
tion of the energies of the noise bursts. Since the envelope function is 
normalized such that this energy e is equal to K 2 , letting y = x and 
e = K 2 yields: 

Q(y) = Pr[e > y] = (e /y) a for y > e„ (5) 
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where e is the energy of the lowest energy noise burst included in the 
distribution. 

The result of the ensuing analysis is the calculation of the conditional 
error rate, in terms of the average number of errors per noise burst, as a 
function of a signal-to-noise ratio. This ratio is defined as the average 
signal energy per data bit divided by the minimum burst energy e . 
For a given modulation system and signal power, a value of Eo may be 
chosen that is sufficiently low to insure that all noise bursts which cause 
errors must have greater energy. Then for any given channel, an impulse 
counter with its threshold set to count impulses with energy greater 
than Eo can find the number of such noise bursts per data bit trans- 
mitted. The product of this number and the conditional error rate is 
then the over-all error rate, in bits in error per bit transmitted. 

The second random variable associated with each noise burst is the 
phase \j/ of the carrier, which is assumed to be uniformly distributed in 
the interval between and 2t. 

The third random variable associated with each noise burst is the time 
of occurrence. Let the sampling instant for an arbitrarily chosen data 
pulse be designated as the time origin t = 0, and let the time of occur- 
rence of the closest noise burst, past or future, be designated as t = t. 
It is assumed that the average interval between noise bursts is so much 
greater than the duration of the bursts, that the probability of two of 
them overlapping is negligibly small. Then, if an error does occur due to 
noise at the sampling instant t = 0, only the nearest burst, occurring at 
t = t, has contributed toward it. For the average rates at which the 
noise bursts have been observed to occur, this assumption is valid. 

For the purposes of the subsequent analyses, the density function 
p(r) is only of interest for absolute values of t less than the duration of 
the noise burst. If the closest noise burst occurs at a time r which is more 
than a noise burst duration removed from the sampling instant, then no 
remnant of the noise is present at the sampling instant, and hence no 
error can occur. This duration is so short compared to the average in- 
terval between bursts that p(r) is essentially constant, equal to p(0), 
over that range of values. Further, for the distributions which have been 
measured, p(0) is approximately equal to the average number of bursts, 
per unit time, /3, of greater energy than e . (This can be shown to be 
exactly true for the case when the intervals are exponentially distributed, 
with a Poisson distribution for the number of noise bursts occurring in a 
given period. It is very nearly true for the other distributions proposed.) 
For absolute values of t greater than the duration of the noise burst, 
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the density function p(r) does not enter the analysis, and no further 
assumptions are required. 

In the following sections, the various modulation systems are ana- 
lyzed, and their performance in the presence of the noise described above 
is determined. 

III. AMPLITUDE MODULATION SYSTEMS 



3.1 Double-Sideband AM with Coherent Detection 

The ideal data receiver using coherent detection consists of a band- 
pass filter H c (co), followed by a coherent detector and a post-detection 
low pass filter H L (u>). When double sideband transmission is utilized, 
the band-pass filter is symmetrical about the carrier frequency a> c . The 
detector multiplies the received signal by a carrier signal with the same 
frequency and phase as that of the modulated carrier at the transmitter. 
The post-detection filter removes those components of the resulting 
signal centered about | o> \ = 2w c , in addition to shaping the resulting 
baseband signal. The post-detection filter is followed by a synchronous 
decision device which samples its input at discrete sampling instants T 
seconds apart. A block diagram of the receiver is shown in Fig. 3. 

At the input to the decision device, the data signal y(t) consists of a 
sequence of identical, except for polarity, pulses y (t) spaced T seconds 
apart 



t(t) = £ a<EVo(t - W) 



(6) 



where a, equals +1 if the ith bit is a mark, and —1 if it is a space. At 
the jth sampling instant t = jT, the data signal y(jT) should equal 
djE with no intersymbol interference present from any other pulse. The 
signal yo(t) is band limited and cannot be restricted to a time slot T 
seconds in width. Each pulse extends over several sampling instants. In 
order to avoid intersymbol interference, yo(t) is constrained to be a 



n(t) 



s c (t) 



*i 2 



H r (w) 



cos((D C t+e) 
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s r (t)+m r (t) s d Ct)+m d (t) y(t)+m(t) 

Fig. 3 — AM receiver with coherent detection. 
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member of a class of functions equal to 1 at t = and equal to zero at all 
other integral multiples of T, satisfying Nyquist's First Criterion. 8 

In the absence of noise, each sample may have one of only two possible 
values. The ith. sample may be ±E, depending on whether the zth bit 
is a mark or a space. When noise is present, the sample may have one of a 
range of values. The decision device is a maximum likelihood estimator, 
producing a mark symbol if the sample value is positive, and a space 
symbol if it is negative. At any sampling instant, the noise will cause an 
error if it has an absolute value greater than E and polarity opposite to 
that of the data signal. The probability of this event occurring, as a func- 
tion of the signal-to-noise ratio, will now be found. 

Let t = denote the sampling instant for an arbitrarily chosen data 
pulse, and let t = r denote the time of occurrence of the closest noise 
burst. As discussed, the probability of two bursts overlapping is negli- 
gibly small. If an error does occur at the sampling instant, it is the result 
only of the closest noise burst. As described in Section II, the noise 
burst n(t) is given by (1), and has a spectrum given by 

N(a) = (K/2){e # iVo(a. - u e ) exp [-j(u - u c )t] 

(7) 
+ e } *N (co + co c ) exp [—j(a + co c )r]}. 

The spectrum of the noise at the output of the post-detection filter 
H L (u) is 

M(w) = (K/4:)N (u>)[e M - e) H c (a) + u c ) 

+ e- il *- n H e ( a - ajy^Hifa). 

As described, the band-pass filter H c (u) is symmetrical about the car- 
rier frequency | co | = w c , therefore, 

H e (u + co c )H L (o>) = H e (u> - co e )# L (co). 

An "equivalent receiving filter" H{w) is now defined as 



so that 



H(a) ^ H e {<* + o>c)H L {u) = H c (a> - u e )H L (<»), (9) 



M(w) = (K/2) cos <t>N (u)H (a) exp (-jut) (10) 



where 4> = \p — 6 is uniformly distributed in the interval between and 
2ir. Let m (t) be defined as the response of the equivalent receiving filter 
H(co) to the normalized envelope no(t), given by 

«o(0 = 7T / No(a)H(u) exp (jut) do. (11) 
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Then, at the output of the post-detection filter, the noise has the spec- 
trum 

M(w) = (K/2) cos<£M„(co) exp (-jut). (12) 

At the sampling instant t = 0, the noise at the input to the decision 
device is equal to 

m(0) = (K/2) cos (j>m (-T). (13) 

An error will occur if the absolute value of w(0) is greater than E and 
if the polarity of the noise is opposite to that of the signal. The probabil- 
ity of error at the sampling instant t = 0, given that the closest noise 
burst has occurred at t = t with relative phase <j>, is thus equal to 

Pr[error | r,0] = \ Pr[m 2 (0) > E 2 | r,0]. (14) 

Substituting the right-hand side of (13) for m(0), a rearrangement of 
terms yields 



r 4E 2 , ' 

Pr[error | T,<f>] = \ Pr e > — 5— — t-. r T,<f> 

1 L cos 2 0m o 2 ( — r) 



(15) 



By (5), this is 

r> r 1 .i 1 feo cos 2 0w o 2 (-t)T , . 

Pr[error | rrf] = - I ^ J . (16) 

The relative phase is uniformly distributed in the interval between 
and 2ir. Therefore, the probability of error at the sampling instant 
t = 0, given that the closest noise burst has occurred at t = t, is given 
by 



„ r 11 C a reoW 2 (— T) 

Pr[error | r] = - [ ^ 



(17) 



C« = ir (^ (cos 2 ^>) a ^. (18) 

l-K JO 



where 

2ir Jo 
The probability of error, at any arbitrary sampling instant, is 

Pr [error] = ^ [^J £ m 2a ( -r)p(r) dr. (19) 

The density function p (r) is only of interest for the range of r less than 
the duration of the noise burst m (t). For larger values of r, m 2a (— t) 
is essentially equal to zero, and the integrand in (19) is zero. As de- 



DIGITAL TRANSMISSION 1709 

scribed in the discussion of the noise model, the density function is es- 
sentially constant and equal to the average number of bursts per unit 
time /3, for the range of t of interest. The probability of error is therefore 
essentially equal to 

Pr[error] = & [^J £ W( -r ) dr. (20) 

The system is constrained to operate with a limitation on the average 
power of the signal on the channel. Under such a limitation, the error 
rate should be given as a function of the average power of the transmitted 
signal, rather than of the peak value E. It may be shown that for the 
double-sideband amplitude modulated signal s c (t), which yields the 
signal y (t ) given by (6 ) , the average power on the channel is 

Fo(co) 2 



w 



- £1 f 



Hip) 



do: (21 



where Yo(oi) is the spectrum of the individual data pulse yo(t). Let 7i 
be defined as the integral 

Fo(«) 2 



r,4 W 

1 2irT lc 



du>. (22) 



loo I ff(a>) 

Then substituting in (21) and (22) into (20), 

Pr[error] = ^ feT f m 2a (-r) dr. (23) 



Substituting t = — r in the definite integral and rearranging some of 
the terms, (23) yields 

p r[err0 r] = ?£f [i /, fj J £ BW(0r f , (24) 

where PFT/eo is the signal-to-noise ratio defined previously as the signal 
energy per bit divided by the energy in the minimum noise burst. The 
rearrangement of terms in the above expression normalizes the integral 
to a dimensionless constant. 

Since F (co) is specified by the system designer, and iV (w) is known, 
the equivalent receiving filter characteristic H (to) may be optimized 
in the sense that the probability of error is minimized. For each value of 
a, however, a different characteristic is optimum. Since the data system 
is to be used over randomly selected channels, with all possible values of 
a, it appears reasonable that it should be optimized in a minimax sense. 
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The system design should therefore be optimized for the worst case, 
which corresponds to a value of a. equal to 1 . If the system then operates 
with an acceptably low error rate over a channel with a value of a equal 
to 1, it will operate with an even lower error rate over another channel 
with a higher value of a. The optimization procedure therefore consists 
of finding the filter characteristic which minimizes the probability of 
error for a value of a equal to one. For that case, the probability of error 
is equal to 



r» r 1 P e « f \ ^o(w) , f 2/.x ,. 

Pr error = — — / ... da> / mo it) at 

where, by Parseval's theorem 

[ mo\t) dt =i [ |Mo( w ) | 2 dco 

J-oo 2.1T J-oo 

f m 2 (t) dt = ±- [ | AT («)ff(«) | 2 do.. 

J— oo .S7T •'-00 



(25) 



(26) 



Then 
Pr[error] = 



32tt 2 LTTtJ ■*-< 



1V«) 



ff(«) 



dw f |Ar (co)^(co) \ 2 du. (27) 

J — CO 



By Schwarz's inequality 



[ F (o))iVo(co) dw [ Si / I ^r f rf " /" I N (a>)H(a>) | 2 d 

J-oo J-oo // («) J-oo 



The equality is satisfied when 

FoU) 
tf( w ) 



= CNo*{u)H*(a:) 



(28) 



(29) 



where C is any real constant and the asterisk denotes the complex con- 
jugate. Since the left hand side of the inequality is independent of 
H(u), it represents the minimum value the right hand side may take. 
The probability of error is therefore minimized when the equality is 
satisfied, and this occurs when the filter characteristic satisfies the relation 



H(o) 



1 Fo(«) 

CN *(o>)' 



(30) 



The phase characteristic of H(u) has no effect on the probability of 
error for the limiting case when a is equal to 1 . The system performance 
is therefore evaluated for zero (or linear) phase shift. The phase char- 
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acteristic does affect the error probability when a is greater than 1 ; this 
effect is considered in detail in a subsequent section on complementary 
delay niters. 

The probability of error is dependent on two parameters of the noise 
bursts, the distribution of their amplitudes-defined by a, and their rela- 
tive frequency — defined by /3. This second dependence consists merely 
of a direct proportionality, and need not be earned along in the computa- 
tion. The performance measure may be normalized by considering the 
conditional error rate, defined as the average number of bit errors per 
noise burst. The probability of error is equal to the average number of 
bit errors per bit transmitted. If this number is divided by the average 
number of noise bursts per bit transmitted, @T, the result is the average 
number of bit errors per noise burst, 

£v _ Pr[error] 



0T 
given by: 

*-tG / *w?J£ bwwi '*- (3i) 

In order to obtain some numerical results, the general expression for 
N given above is to be evaluated for a special case which is of interest. 
In the transmission of data over narrow-band telephone channels, the 
spectrum of the individual data pulses y (t) is often the "raised cosine" 
spectrum given by 

Fo(«) = (T/2)[l + cos (uT/2)] for | « | < 2t/T. (32) 

This signal satisfies all three of Nyquist's criteria: 



(l.) yodT) 



1 for i = 
Ofort = 1,2, 3, 



(3.) The envelope of y (t) approaches zero very rapidly as | t | increases. 

For the sake of simplicity it will be assumed that in the narrow band of 
interest | w | < 2ir/T, the spectrum of the noise No(u) is constant. Since 
the energy of the noise burst n (t) in the band of interest is equal to 2 
watt-seconds, the noise spectrum is given by 

N (w) = y/f for | « | < 2r/T. (33) 
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The optimum filter characteristic is given by 

H(u) = {Ml + cos (a>7y2)]}* for |»| < 2t/T 



ff(w) = cos (ojT/4) 



for | to | < 2ir/T. 



(34) 



Note that the optimization permitted an arbitrary scale factor for H(u), 
since signal and noise would be identically affected. This scale factor 
has been chosen such that the maximum value of H(co) is unity. 

The expression for N has been numerically evaluated on a digital 
computer, for values of a equal to 1, 2, and 3, chosen to be representa- 
tive. The results are tabulated here, and are plotted in Fig. 4. 
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Fig. 4 — Performance of double-sideband AM system with coherent detection. 
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3.2 Single- Sideband AM with Coherent Detection 

The ideal data receiver has the same form as the one described in the 
previous section. The only difference is that when single-sideband trans- 
mission is utilized, the band-pass filter H c (u) does not pass frequencies 
below co c 

H e (u) = for !«!<«,. (35) 

At frequencies above u) c , the band-pass filter is identical to that for the 
double-sideband case. 

As in the previous analysis, let t = denote the sampling instant for 
an arbitrarily chosen data pulse, and let t = r denote the time of occur- 
rence of the closest noise burst. As in the double-sideband case, the 
spectrum of the noise burst at the output of the post-detection filter is 
given by (8). However, since H c (o>) is equal to zero for all | co | less than 
Wc , this spectrum is equal to 

M(u) = (K/4)M (co) [cos0 + jsgn(u) sin 0] exp (—Jut), (36) 



and the noise burst is 

m(t) = (K/4:)[m (t - t) cos<f> — n\ (t — t) sm<f>] 
where mo(t) is the Hilbert transform of m (t), given by: 



™o(t) = ~ ji- sgn(u>)N (u>)H(a>) exp (jut) day. (37) 

2tt J-oo j 

At the sampling instant t = 0, the noise at the input to the decision 
device is equal to 

m(0) = (K/4)[m (- T ) cos - iho(-r) sm<f>]. (38) 

Substituting the right-hand side of (38) for w(0) in the expression for 
the probability of error, (14), and rearranging terms, 

Pr[error | t,4>] 

16E 2 



- *Pr<e > 



M-t) cos — lho(— t) sin 4>f 
By (5), this is equal to 



T,<f> 



(39) 



-r, r I ,-, 1 fe [rao(— t) cos </> — ra (— t) sin <f>f\ 
Pr[error | t,0] = - < ^q^ f » ' 40 ' 

and may be rewritten in the form 
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Pr[error | r,«] = \ [j^T [m 2 (-r) + m,\-r)] a (cos 2 x )" (41) 

where, for any given value of r, the angle 

_im (— t) , . v 

X = + tan — ? - (42) 

wo(-r) 

is uniformly distributed in the interval between and 2x. Averaging 
(41 ) with respect to x and r yields 

Pr[error] = ^= [^ J £ [W( -r) + *fto 2 (-r)] a dr. (43) 

It may be shown that, for the single-sideband amplitude-modulated 
signal s c (t) which yields the signal y(t) given by (6), the average power 
on the channel is equal to 

W = ±E% (44) 

where 7i is the integral denned by (22). Substituting (44) in (43) and 
rearranging some of the terms yields: 

Pr[ error] =^f[\h ^J £ PW(0 + Tm 2 (t)] a % (45) 

where WT/eo is the signal-to-noise ratio defined previously. 

In the same manner as for the double-sideband system, the equivalent 
receiving filter characteristic H{o>) is designed to minimize the prob- 
ability of error for the case when a is equal to 1. For that value of a, the 
error probability is given by 



(46) 



Pr[error] 

- L[wi] /1|f$H/-> !(( > * + £><» 4 

By Parseval's theorem 

I m 2 (t)dt = f m Q \t)dt = i- [ \N (u)H(o>)\ 2 du>. (47) 

J — oo -J— oo 41T J— co 

Substituting (47) into (46) yields the same expression as for the 
double-sideband case, given by (27). The optimum filter characteristic 
is therefore the same, satisfying (30). 

The conditional error rate, defined as the average number of bit errors 
per noise burst, is given by 
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L«ri 7 -^-Tf 

2 [_4 l pptJ J- 



\Tmf(t) + Mf{t)r*. (48) 



In order to obtain numerical results, the general expression for N is 
evaluated for the special case described previously. The spectrum of an 
individual data pulse is the raised cosine spectrum, and the noise spec- 
trum is assumed to be constant in the band of interest. The expression 
for N has been numerically evaluated for values of a equal to 1, 2 and 3. 
The results are tabulated here, and are plotted in Fig. 5. 

N 



0.124 (eo/WT) 
0.0241 (e /WTf 
0.0069 (eo/WT)* 



3.3 AM with Envelope Detection 

The ideal data receiver using envelope detection consists of a receiving 
filter H c (u), symmetrical about the carrier frequency o c , followed by 



m iQ-2 



u. |0" 4 



2 

10-5 

5 



a = \ 

,2 



^ = 10 LOG 10 [^p-\ IN DECIBELS 



Fig. 5 — Performance of single-sideband AM system with coherent detection. 
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an ideal envelope detector. The output of the envelope detector is applied 
to a synchronous decision device. A block diagram of the receiver is 
shown in Fig. 6. The signal y{t) at the input to the decision device has 
the form of (6 ) . At the ith sampling instant, it is equal to 

yiiT) = |a,|. 

Since polarity information is lost, and only magnitude information is 
delivered to the decision device, unipolar signaling must be used. If the 
ith bit is a mark, a» is equal to +1; if the ith bit is a space, a t is equal 
to zero. The threshold of the decision logic is set at a level \E, where X 
has a value between zero and one. If the ith sample is greater than XE, 
a mark symbol is produced; if less than XE, a space symbol is produced. 
The value of X is chosen to minimize the probability of error in the 
presence of noise. 

Let t = denote the sampling instant for an arbitrarily chosen data 
pulse, and let t = r denote the time of occurrence of the closest noise 
burst. The noise burst on the channel is given by (1). At the input to 
the detector, the noise burst is given by 

■ r (t) = Km (t - r)[cos<£ cos (u c t + 0) - sin sin (<a t + 0)] (49) 



m 



where the previously described phase difference 

= ^-0, 

is uniformly distributed in the interval between and 2ir. 

At the input to the detector, the combined signal and noise voltage 
is equal to 

Vi(t) = [y r (t) + Km (t - t) cos^] cos (a> c t + 0) 



— Km (t — t) sin sin (io c t + 0). 



(50) 



The detector output is 



Vo(t) = V[y r (t) + Km (t - r) cos <j>f + [Km (t - r) sin <j>f. (51) 



n(t) 



s c (t) 



H>) 



ENVELOPE 
DETECTOR 



DECISION 
DEVICE 



v l (t)=s r (t)+m r (t) v d (t) = y(t) + m(t) 

Fig. 6 — AM receiver with envelope detection. 
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If, at the sampling instant t = 0, a space has been sent, the voltage 
at the input to the decision device is equal to the noise alone 

v (0) = \Km (- T )\. (52) 

An error will occur if this magnitude is greater than \E. The probability 
of error at the sampling instant t = 0, given that a space has been sent 
and that the closest noise burst has occurred at t = r, is therefore equal 
to 



Pr[error | space, r] = Pr[erao 2 (— r) > \ 2 E 2 
and, by (5), this is 

eom ( — i 



(53) 



-n r i i reom 2 (— t)T 

Pr[error | space, t] = I — — — I 



(54) 



If, at the sampling instant t = 0, a mark has been sent, the voltage 
at the input to the decision device is equal to 



vo(0) = +V[E + Kmoi-r) cos 0] 2 + [Km (- T ) sin <f>f (55) 

and an error will occur if this is less than \E. This is shown graphically 
in the phasor diagrams of Fig. 7. An error will occur if the vector sum of 
the signal E and the noise Km ( — t) making an angle <f> with E falls 
within the circle of radius \E. 

Fig. 7(a) shows the phasor diagram when ra (— t) is positive. An 
error can occur only for values of <£ in the interval 

■k + sin -1 X ^ (f> ^ 7r — sin -1 X. 

For all other values of <f> the vector sum cannot fall within the circle. 
When <f> is within the above interval, an error will occur if Km (—T) 




Fig. 7 — (a) Phasor diagram when m (—r) is positive, (b) Phasor diagram when 
?«o(— t) is negative. 
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lies in the interval 

m 2 > Km (—T) > mi (56) 

where wii and W2 are functions of 0. By the Law of Cosines, m x and mi 
must satisfy 

X 2 ^ 2 = E' 2 + m - 2mE cos (ir - 0), (57) 

hence, 

(1 - X 2 )£ 

cos (t — 0) + VX 2 — sin 2 (it — 0) 

(i - V)E (58) 



mi = 



W-2 



COS 



(tt - <f>) - VX 2 - sin 2 (t - 0) 



The probability of error at the sampling instant t = 0, given that a 
mark has been sent and that the closest noise burst has occurred at 
t = t, such that rtioi—T) > 0, with phase difference <j>, is equal to the 
probability that the noise is within the range (56), and is 

Pr[error | mark, t,v: (— t) > 0,0] 

[_m 2 ( — t) m \ — T)\ J 

for x + sin - X > <p > r — sin - X 

and is equal to zero for all other values of 0. By (5 ) this is equal to 

Prferror | mark, t, ra (— t) > 0,0] 

'e m 2 (-r)[cos (tt - 0) + VX 2 - sin 2 (7r - 0)] 2 1 



(1 _ X 2)2.E 2 

e m 2 (— t) [cos (it — 0) — \/\ 2 — sin 2 (t — 0)] 2 V* 
(1 - X 2 ) 2 £ 2 J 

for 7r + sin - X > > t — sin - X. 

The phase difference <f> is uniformly distributed in the interval between 
and 2tt. Therefore, averaging (60) with respect to <f> yields: 

Pr[error | mark, t, ra () (— t) > 0] 

— [cos — \/\ 2 — sin 2 0] 2 " J d<f>. 
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Fig. 7(b) shows the phasor diagram when m (— t) is negative. By 
reasoning similar to that presented above, it can be shown that the 
probability of error at the sampling instant t = 0, given that a mark has 
been sent and that the closest noise burst has occurred at t = t, such 
that m (— t) < 0, is also given by (61). The probability of error 
given that ra (— t) is positive is identically equal to the probability of 
error given that it is negative. Since the two events are disjoint 

Pr[error | mark, t] 

i r ^ ( \ ~\ a /* H ' ,|— '^ 
= Alf^m\l lteos* + Vx^li^r (62) 

— [cos — \/x 2 — sin 2 <f}f a } d<f>. 

Marks and spaces are assumed to occur with equal probability. There- 
fore, the probability of error at the sampling instant I = 0, given that 
the closest noise burst has occurred at t = t, is equal to one-half the 
sum of (54) and (62). Let a signal-to-noise ratio coefficient be defined 
by 



(63) 



/ " (X)= {^ + ,(l- 1 X^ f" 1X ^ 0S * + Vx 2 -sin 2 rf" 

- |cos <f> - Vx 2 - sin 2 0p} d<t>\ l ' a . 

Then, the probability of error, given t, equals 

df ii 1 pa(X)r m 2 (-r)T 

Pr[error \ t\ = - — . (64) 

As shown previously, averaging with respect to t yields that the prob- 
ability of error at any arbitrary sampling instant is essentially equal to: 

Pr[error] = \ [^^J £ W( -f ) dr. (65) 

The value of the threshold level X is to be selected to minimize the coef- 
ficient /a(X). The functions / a (X) and df a (\)/d\ have been evaluated 
for values of a equal to 1, 2, and 3. These are plotted in Figs. 8 and 9. 
respectively. Since a value of a equal to 1 yields the highest error rate, 
the threshold level X should be selected to minimize /i(X). The error 
rates for other values of a will always be lower than that for a value of 
a equal to 1. However, /i (X) exhibits the broadest minimum and is most 
amenable to compromise. As may be seen in Fig. 8, a threshold level X 
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Fig. 8 — Signal-to-noise ratio coefficient. 
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Fig. 9 — Rate of change of signal-to-noise ratio coefficient. 
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equal to 0.6 is a satisfactory compromise. At that value of X, fa equals 
3.65, fa equals 3.55, and fa equals 3.65. 

It may be shown that for the unipolarly keyed modulated signal s c (t), 
the average power on the channel is equal to 

w = e 2 r~/i + — f ° (o) n . (66) 

L8 8T 2 tf(O) I J 

The second term is recognized as the result of a carrier frequency com- 
ponent which is present because the average value of the baseband data 
signal is not equal to zero. Substituting (66) into (65) and rearranging 
some of the terms, 



Prferror] = f Pl 



7 1 + 77T 2 



yo(0) I I i r '- ,/, 

HW\] \ ]_ [Tm 2 (t)) a ^ (67) 



8WT 



where WT/bq is the signal-to-noise ratio denned previously. 

In the same manner as for the previously considered systems, the filter 
characteristic H(u) is selected to minimize the error probability when a 
is equal to one. For that value of a, the probability of error is given by 

Pr[error| = ^rdLiw 1 + t 8(w t w / x 

(68) 
• [ \NMH(a) | 2 do,. 

J— as 

By Schwarz's inequality, the above expression is minimized when the 
amplitude characteristic of H((a) satisfies the relation 



HM |2 Yo(u) 



No*(«) 



[l + ^ *(«)]'. (69) 



The requirement of an impulse at zero frequency is equivalent to carrier 
suppression at the transmitter. The infinite gain at the receiver then 
restores the dc component of the base-band data signal. Since the 
construction of such a filter is not feasible, it is assumed that the filter 
H(u) has the suboptimum "smooth" characteristic given by (30). 

The conditional error rate, defined as the average number of bit errors 
per noise burst, is 



In order to obtain numerical results, the general expression for N is 



»-*£'• 



A + * 
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evaluated for the special case described previously. The spectrum of an 
individual data pulse is the raised cosine spectrum, and the spectrum of 
the noise is assumed to be constant in the band of interest. The expression 
for N has been numerically evaluated on a digital computer, for values 
of a equal to 1, 2 and 3. The results are tabulated below, and are plotted 
in Fig. 10. 

N 



0.455 (e /W T) 
0.456 (e /WTf 
0.586 (e /WT) 3 



IV. FREQUENCY SHIFT KEYING SYSTEM 



The ideal data receiver for a frequency shift keying system consists 
of a receiving filter H c (ui) which is symmetrical about the carrier fre- 
quency Qi c , followed by an ideal frequency discriminator. The output of 
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Fig. 10 — Performance of AM system with envelope detection. 
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the discriminator is the instantaneous frequency w,- of the signal (plus 
noise) at its input. The discriminator is followed by a synchronous de- 
cision device identical to the one in the coherent AM systems. A block 
diagram of the receiver is shown in Fig. 11. 

Sunde 9 has shown that in the absence of noise, intersymbol interfer- 
ence can be eliminated even though the signal at the input to the de- 
tector must be bandlimited. The frequency shift signal is generated by 
mixing the outputs of two synchronized oscillators; one oscillating at a 
mark frequency co m , the other at a space frequency co s . For this discus- 
sion, let the mark frequency be higher than the space frequency; the 
choice is arbitrary. The carrier frequency is defined as the mid-frequency 

C0 C = (w„, + W s )/2 

and the shift frequency is defined by : 

2o)j = o) m — a)j . 



The modulated signal is 

s(t) = i{.[l-+ git)] cos [(«« + m)t + 



- [1 ~ g(t)][C08 (co c - a> d )t+ 0]. 
At the sampling instants, the mixing function is equal to 



(71) 



giiT) = ±1 

so that the instantaneous frequency at the sampling instants is 

m(iT) = Wc ± &d 

depending on whether the ith bit is a mark or a space. 
The modulated signal may be rewritten as 

s(t) = sin w d t sin (o3 c t + 9) — g(t) cos wd cos (<a c t + 6). (72) 

After successive filtering, amplification, and transmission, the signal at 
the input to the discriminator, in the absence of noise, is equal to 



s c Ct) 



H c (cu) 



s r (t)+m r (t) 



DISCRIM- 
INATOR 



U,(t) 



DECISION 
DEVICE 



Fig. 11 — Receiver for frequency shift keyed system. 
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8 r (t) = P(t) cos (a, c t + 0) - Q(t) sin (aj + 0), (73) 

where 

F0) - -, W 
Q(<) = — -Ssin w d tf. 

As a further restriction, the shift frequency is made equal to the bit 
rate, so that 

m = t/T 

and 

y(t) = ^(-lYaiEyoit-iT), (75) 

i 

where, as before, a, is + 1 if the ith. bit is a mark and — 1 if it is a space, 
and yo(t) satisfies Nyquist's First Criterion. Sunde has shown that these 
conditions yield a signal which in the absence of noise presents no 
intersymbol interference. 

As in the analysis of the AM systems, let t = denote an arbitrarily 
selected data sampling instant, and let t = r denote the time of occur- 
rence of the closest noise burst. The noise at the output of the receiving 
filter is given by (49). With both signal and noise present at the dis- 
criminator input, the output is 

where 

x(t) = Kmo(t — t) cos<f> 

y(t) = Km (t — t) sia<f) 
and is equal to 

_ [Pit) + x(t))[Q(t) + y(Q] - [Q(t) + y(t)][P(t) + x(t)] (78) 
[P(0 + x(t)f + [Q(t) + y{t)f 

Since the denominator of co,-(0 is always positive, the decision device 
will produce a mark if 

V = [-Ea Q 4- jKwio(-t) cos0][- (Ett/T) + Km a {-r) sintf>] 

(79) 
- [Km (-T) sm<t>][P(0) + Krho(— t) cos<£] 

is positive, and a space if it is negative. 



(77) 
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For the cases of interest, it can be shown that P(0) is considerably 
less than x(0) for those noise bursts which cause errors. In addition, 
P(0) is equal to zero. For these reasons, the P(0) term is dropped in 
the subsequent steps, since to retain it would unnecessarily complicate 
the analysis. 

When a mark has been sent, and a is equal to +1, V is negative and 
an error occurs when 



K[(T/T)rh (-T) sin<£ + w (-r) cos 4>] > E. 



(80) 



The probability of error at the sampling instant t = 0, given that a 
mark has been sent and that the closest noise burst has occurred at 
t = t with phase difference <f>, is equal to : 



Pr[error | mark,r,0] 



-i* 



E 2 



e > F 



rp 

— m ( — t) sin + m Q ( — T) cos <f>] 



T,0| 



(81) 



and by (5), this is equal to 
Pr[error | mark,T,0] 



--= -< 
2 



eol - Wo( — t) sin 4> + mo( — t) cos <f> 



II 



E- 



(82) 



When a space has been sent, and a is equal to —1, V is positive and 
an error occurs when : 



K[(T/ir)m Q (-T) sin - m (-T) cos 0] > E. 



(83) 



Following the reasoning presented above, the probability of error at the 
sampling instant t = 0, given that a space has been sent and that the 
closest noise burst has occurred at t = r with phase difference <f>, is 
equal to 



Pr[error | space,r,<£] 



Bo — Wo( — t) sin <f> — m ( — T) cos <f> 
It 



I) 



E"- 



(84) 



Marks and spaces occur with equal probability. Therefore, the prob- 
ability of error, given r and <f>, is 
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Pr[error | r,0] = i[j|J |[^ m (-r)J + [moC-r)] 2 !" 



cosU+ tan 7 — 

1 L t m { — 



ar 



(85) 



-f- scos' 



[— <";&$$)■ 



Averaging with respect to 4> and t yields 



Pr[error] = -^ \ ^- 



=[!]'£§•«->] 



-mo(-r) + [m (-r)ry dr. (86) 



It can be shown that, for the frequency shift keyed signal s c (t), the 
average power on the channel is 

W = E 2 [(l/±Ao 2 ) + tfj, (87) 

where 

A = |ff(±«.)| = \H(±a m )\. (88) 

The first term in the bracket of (87) is the result of discrete components 
of the signal at the mark and space frequencies. Substituting (87) into 
(86) and rearranging some of the terms yields 



Pr[error] = 



0C a T 



eo 



2WT J 



(89) 



£|r[^i* (oj+r[m.«)] , r 



where W T/e is the signal-to-noise ratio defined previously. 

In the same manner as for the systems considered previously, an 
equivalent receiving filter characteristic H(w) is to be found which 
minimizes the error rate when a is equal to 1. For that value of a, the 
probability of error is equal to 

Pr[ error] = &. \ «L~\ \ J- + ±- ( " I !• W 
iiienoij g [^ r J[ 2 Ao 2 ^ 2kTI*\H 



2 "I 

do> 



H(«) 
j^ J~W(0 +m 2 (f) eft, 



(90) 



where 



j- = i r r 

A 2 2 1=0 L 



5 ( co - ^ ) + 8 



+ - 



ff(co) 



eta. 



(91) 



(92) 
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By Parseval's theorem 

/ rho 2 (t) dt = — / w 2 | MM | 2 du. 

J— oo 27T J— oo 

Therefore, for a value of a equal to 1, (90) may be rewritten 

+ 1 ^o(«> ri I tf^t r *•- r fs « f + 1) i ^o(«)^(«) i 2 d». 

J | //(,0J_) I •'-oo \7T" / 

For ease of notation, let 

o(«) ^ Tr/2{«[«- (tt/7 , )]+ 5[co+ (tt/T)]) + | Fo(co) | 8 

&( w ) ^ [(T7»V + 1] | Ar (co) | 2 (93) 

r(w) = |//(w) | 2 . 

To minimize the probability of error, it is necessary to find the function 
h((j)) which minimizes the product. 

P = [*?!&!>&,(* b(fi»)h(. u ) dco. (94) 

.Loo AI(CO) J-oo 

The minimum occurs when the variation 
8 p = ( ^M do; ( b(u)Sh(o») du 

J-00 At(ti)) J-oo 

_ f Z^shM du> [ b(o)h(u) du> (95) 

J-ao ft Z (a>) JLoo 

8P = f f [^4 b M - ^ &(«l)AWlaW d«i tfu> 2 

J-=oJ-» L^(wi) « 2 (w«) J 

is equal to zero. Since it must be equal to zero for any variation 8h(<at), 
the bracketed term must be zero 

^ h W - gH & ("0 * to = (96) 

for all values of u»i and a>2 . This requires that h 2 (a))b(a))/a(w) be equal 
to a constant. It can be shown that the second variation 5 2 / J is positive 
when the above condition is met, so that P is truly at a minimum, and 
not at a maximum or stationary point. To minimize the probability of 
error, the filter characteristic must satisfy 
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H(w) 



r ^[«(«-^) + *(« + ^)]+|r.(«) 



s— ) 



iVo(co) 



(97) 



The impulses in the receiving filter characteristics correspond to sup- 
pression of the discrete components at the transmitting filter. Since it 
is not feasible to construct such a filter, it is assumed that the sub- 
optimum "smooth" filter, given by 

I Y (u) 1 
I H(») | 2 



g.- + .y 



N (a>) 



(98) 



is to be used. 

The conditional error rate, defined as the average number of bit errors 
per noise burst, is given by 



» s-r*-T ri- 



ff = 



2 IWTj \_±A 



+ o 



1 ~\ a 

2 7l J 

•r 

J— 00 



[T "I 2 'i « -// 

-fftoCOj +T[m(t)] 2 } ,„. 



T - 



(99) 



r 



The general expression for iV is evaluated for the same special case 
described previously. The spectrum of an individual data pulse is the 
raised cosine spectrum, and the spectrum of the noise is assumed to be 
constant in the band of interest. The expression for N has been numeri- 
cally evaluated, on a digital computer, for values of a equal to 1, 2, and 
3. The results are tabulated below, and are plotted in Fig. 12. 

N 



0.402 eo/WT 
0.392 (e<,/WT) 2 
0.48 (e /WT) 3 



V. PHASE SHIFT KEYING SYSTEM WITH DIFFERENTIALLY COHERENT DE- 
TECTION 



5.1 Binary (Two Phase) System 

The ideal data receiver utilizing differentially coherent detection con- 
sists of a receiving filter H c (co) symmetrical about the carrier frequency 
oic , followed by a detector. The differentially coherent detector multi- 
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Fig. 12 — Performance of frequency shift keyed system. 

plies the received signal by the signal which had been received one bit 
duration earlier. The output of the detector is applied to a synchronous 
decision device of the type described previously, in the analysis of the 
double-sideband AM system. A block diagram of the receiver is shown 
below, in Fig. 13. 

In the absence of noise, the signal at the input to the detector is of 
the form 

s r (t) = y(t) cosojc* (100) 

where the modulating signal y(t) satisfies (6) The multiplier a, may be 
±1. If the ith bit is a mark, a, is made equal to ai_i ; if a space, a, 
equal to — a,_i . In the absence of noise, the signal at the output of the 
detector is equal to 

v (t) = y(t)y(t - T) cos u e (t) cos u e (t - T). (101) 

The carrier frequency w c is selected to be an integral multiple of the bit 
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n(t) 




v.(t)=s r (t)+m r (t) 



Fig. 13 — Receiver for binary phase shift keyed system with differentially 
coherent detection. 



rate, so that cc c T is an integral multiple of 2ir. At the ?*th sampling in- 
stant, the signal at the output of the detector is equal to 

v (iT) = y(iT)y[(i - 1)T] = a^^E 2 . (102) 

The sample is equal to -\-E 2 if the ith bit is a mark, and —E 2 if it is a 
space. The decision device produces a mark symbol if the sample is 
positive and a space symbol if it is negative. 

Let t = denote the sampling instant for an arbitrarily chosen data 
pulse, and let t = r denote the time of occurrence of the closest noise 
burst. At the sampling instant t = 0, the output of the detector is 

Wo(0) = [aoE + Kmoi-r) cos 4i\[(Lj8 + Km (- T -T) cos *]. (103) 

To find the probability of error, given t and \f/, it is necessary to find the 
ranges of values of K which cause the polarity of i>o(0) to be reversed. 
For ease of notation, let the following two functions of r and \f/ be 
defined: 



Bi(rj) ^ 



Bt(r,4,) & 



E 



TOo(-t) cos \p 



E 



m (—T — T) cos >p 



The output of the detector, at the sampling instant t 
written 

*(0) = E 2 [ao + (K/B^a^ + (K/Bt)]. 



(104) 

0, may be 

(105) 



Depending on the values of r and \f/, Bi and B 2 may each be either posi- 
tive or negative, so that there are four possible combinations of polarity. 
In addition, either B x or B 2 may have the larger absolute value. Each of 
these eight possible combinations of polarity and relative absolute value 
must be investigated separately. The range of K which causes a reversal 
of polarity, depending on the values of a and a_ x , is to be found. This 
is most clearly done graphically. In Fig. 14, the eight possible combina- 
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tions are shown separately and, for each of the four equally probable 
combinations of a u and a_i , the range of K which causes an error is 
shown shaded. Because the four combinations of o and o_i are equally 
probable, each with probability of occurrence equal to j, the probability 
of error, given t and xf/, for any one of the eight possibilities for B x and 
Z? 2 , is equal to 

Pr[error | r,\f/\ = | £3 Pr[K in shaded region]. (106) 

"o.a-1 

It may be seen in Fig. 14 that, for each of the eight possibilities, this is 
equal to 

Pr[error | t,*} = \ Yx\K > min (\B 1 \,\B i \)\ tM (107) 

The probability of error, given r and \f/, may be rewritten, 

Pr[error | r,\f/] 

E 1 



= h Pr { e > 



r,yj/ 



(108) 



\ cos 2 \f/ max [m 2 ( — r),W( — r — T)] 

and, by (5) this is equal to 

r - , 1 fe cos 2 ^ max [m 2 (— T),m 2 (— t - jT)]\ a , . 
Pr[error | r,$ = ^ \ ^i 1 • ( 109 ) 

Averaging with respect to \p and r yields 

Pr[error] = ^ (V]" £ {max [m 2 (-r),W(-r - 7 1 )])" rfr. (110) 

The signal on the channel s c (t) is the same as that for the double-side- 
band AM system. The average transmitted power is therefore given by 
(21). Substituting this into (110), dividing by /37\ and rearranging 
terms, yields that the conditional error rate is equal to 



N = c f\: 



[i h fjj £ {max [Tm 2 (t),Tm 2 (t - T)}}° * (111) 



where WT/vq is the signal-to-noise ratio as defined previously. 

In order to obtain numerical results, the general expression for N 
given above is evaluated for the same special case as the preceding sys- 
tems. The spectrum of a single data pulse is the raised cosine spectrum, 
and the spectrum of the noise burst is assumed to be constant in the 
band of interest. The equivalent receiving filter characteristic H(u) is 
the same as that for the AM systems. The expression for N has been 
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numerically evaluated, on a digital computer, for values of a equal to 
1, 2, and 3. The results are tabulated below, and are plotted in Fig. 15. 

N 



0.236 (e /WT) 
0.109 (eo/WT) 2 
0.0625 (e /WT) 3 



5.2 Quaternary (Four Phase) System 

A quaternary phase shift keying system transmits a signal which, at 
the discrete sampling instants, may have any one of four possible phases 
spaced 90° apart. Such a signal is mathematically equivalent to two 
binary signals in quadrature, where each binary signal is of the form 
described in the preceding section. 
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Fig. 15 — Performance of binary phase shift keyed system with differentially 
coherent detection. 
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The transmitted signal consists of an "in-phase" binary signal and 
"quadrature" binary signal. The ideal data receiver consists of a receiv- 
ing filter H e (a), centered about the carrier frequency a> c , followed by two 
separate binary receivers. One of these, identical to the receiver of the 
preceding section, is sensitive to the "in-phase" binary signal. The 
other binary receiver, preceded by a phase shifting network with a phase 
characteristic which is equal to — tt/2 radians throughout the frequency 
band of the received signal, is sensitive to the "quadrature" binary sig- 
nal. A block diagram of the receiver is shown in Fig. 16. The system 
operates in the following manner. In the absence of noise, the signal at 
the output of the receiving filter is of the form 

s r (t) = y a (t) cos aid — y b (t) sin w e t, (112) 

where the modulating signals are each of the form given by (6). At the 
sampling instants the second term of (112) is equal to zero, and the "a" 
system is identical to the binary system described in the preceding sec- 
tion. In the absence of noise, the signal at the input to the "6" detector 
is given by the Hilbert transform 

3(0 = Vb cos Wet + y a sin a> c t. (113) 

At the sampling instants the second term of (113) is equal to zero, 
and the "6" system is also identical to the binary system. At any 
arbitrarily chosen sampling instant Pr["a" bit in error] and Pr["6" bit 
in error] are both given by (110). The average number of bits in error 
per symbol transmitted is equal to 

P = Pr["a" bit in error, "6" bit correct] 

+ Pr["6" bit in error, "a" bit correct] (114) 

+ 2 Pr["a" bit and "6" bit both in error]. 

The Venn diagram in Fig. 17 shows that this is equal to 

P = Pr["a" bit in error] + Pr["6" bit in error]. (115) 

The average number of bits in error per symbol transmitted is therefore, 

P = PC |"^T jf Imax [m 2 (-r) ,m 2 (-r - T))\ a dr. (116) 

Since there are, on the average, /37 1 noise bursts per symbol transmitted, 
the average number of bits in error per noise burst is equal to 

N = P//3T. (117) 
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Fig. 16 — Receiver for quaternary phase shift keyed system with differentially 
coherent detection. 

The transmitted signal consists of two binary signals in quadrature. 
Since the two binary signals are orthogonal, the average power of the 
transmitted signal is equal to the sum of the average powers of the two 
binary signals, and is twice the power for the binary system. Substitut- 
ing this and (117), into (116) and rearranging terms yields 



N = C t 



\_ h w\\ L 



[max [Tm 2 (t) ,Tmo (t - T)\ 



^, (118) 



where in this case, WT '/2e is the signal-to-noise ratio, described pre- 
viously as the signal energy per bit divided by the minimum energy per 
noise burst. When N is plotted versus 10 logm [WT/2e ] db as the ab- 
scissa, the value is twice that for the binary case. 




Fig. 17 — Venn diagram for decision outcomes. 
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VI. COMPARISON OF THE MODULATION SYSTEMS 

Iii the preceding three sections, the various data transmission systems 
have been analyzed, and the performance of each in the presence of 
impulsive noise has been determined. It is of interest to compare these 
results and to rank the systems as to performance. For each of the 
modulation systems, an expression has been derived for the average 
number of bit errors per noise burst, as a function of the signal-to-noise 
ratio. The systems may be ranked by comparing the signal-to-noise 
ratios required by the different systems for the same error rate. 

Such a comparison is done here for the special case for which N has 
been evaluated. The spectrum of an individual data pulse is the raised 
cosine spectrum, and the spectrum of the noise burst is assumed to be 
constant in the band of interest. The equivalent receiving filter for each 
system is the one which minimizes N when a is equal to 1. (In those cases 
for which the optimum filter has impulses in its response, corresponding 
to carrier suppression at the transmitter, the suboptimum "smooth" 
filter is assumed. This is described in the preceding sections.) 

For each of the modulation systems, the average number of bit errors 
per noise burst is given by the expression 



*-N? 



(119) 



where the constant C depends on the type of modulation system, the 
characteristics of the receiving filter, and the value of a. The systems 
may be ranked by comparing the values of C. The ratio between the 
values of C for any two systems is the difference in signal-to-noise ratio 
required for equal error rate. In Fig. 18, the values of C for the various 
modulation systems are plotted as functions of a. 

The comparison of Fig. 18 shows the modulation systems to be ranked 
in the following order : 

(1.) Single-sideband AM with coherent detection. 

(2.) Double-sideband AM with coherent detection. 

(3.) Phase shift keying with differentially coherent detection. 

(4.) Frequency shift keying. 

(5.) AM with envelope detection. 
It is interesting to note that this is the same ranking as has been de- 
termined for performance in the presence of Gaussian noise. 

VII. COMPLEMENTARY DELAY FILTERS 

In the sections on the performances of the various modulation sys- 
tems, the general expressions for the error rate have been evaluated for 
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Fig. 18 — Comparison of modulation systems. 



a special case which is of wide interest. For the purpose of those evalua- 
tions, the equivalent receiving filter //(a>) has been given an amplitude 
characteristic A(u) which minimizes the error rate for the worst case, 
when a is equal to 1, and a phase characteristic <j>(<a) which is zero or a 
linear function of frequency 



0(a)) = —Deo. 



(120) 



Such a phase characteristic causes a pure delay of D seconds, with no 
distortion. Since this corresponds merely to shift of the time axis, and 
affects signal and noise identically, it has no effect. It has been pointed 
out previously that, for values of a greater than 1, a phase characteristic 
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0(w) which is a function of frequency other than linear can serve to re- 
duce the error rate. This effect is now discussed. 

When the phase characteristic is other than a linear function, the 
envelope delay 

D( u ) = - f «(«) (121) 

is not constant, but varies with frequency. When the noise burst passes 
through such a filter, its various frequency components are delayed by 
different amounts, and its energy is spread out over many bits. The peak 
value of the spread out burst is much lower than the peak value of the 
original burst, so that many bursts which would have caused errors no 
longer do so. 

In the analyses of the systems, the individual data pulses at the re- 
ceiving filter output are each constrained to be a specific function ijo(t). 
This implies that at the transmitter the data pulses pass through a 
transmitting delay filter, with envelope delay equal to 

D T (u) = D - D(a) 

in the frequency band of the data signal. The constant D, at least as 
large as the maximum value of D{o>), is necessitated by the fact that, 
for realizability, D T (u>) must be positive. The two delay filters exactly 
complement one another, so that the net effect on the data signal is a 
pure delay, with no distortion. The noise burst, occurring on the channel 
after the transmitting filter, only passes through the second filter and is 
spread out. 

The use of such complementary delay filters, based on heuristic reason- 
ing of the type given above, has been suggested previously. 4 ' 11 In this 
section, the improvement which results from the use of such filters, in 
terms of the equivalent increase in signal-to-noise ratio which would be 
required for an equal reduction in error rate, is presented. 

Three types of delay networks, which lend themselves to synthesis, 
have been considered. One of these, called a "linear delay network," 
has an envelope delay characteristic given by 

D x («) = (D m T/2r)oi sgn(aj) for | a \ < 2r/T (122) 

and hence a phase characteristic equal to 

0i(«) = -(Z) m 7y4ir)<o 2 sgn( w ) for | co | < 2tt/T, (123) 

where D m is the maximum delay in the band of interest. The other two 
types are called "sinusoidal delay networks." One of these, with a half 
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cycle of sinusoid in the band of interest has an envelope delay character- 
istic 

D 2 ( w ) = D m cos cor/4 for | <o | < 2r/T (124) 

and hence a phase characteristic equal to 

02 (a,) = _ (±D m /T) sin a.T/4 for | a. | < 2tt/T. (125) 

The third type is a sinusoidal delay network with a full cycle of sinusoid 
in the band of interest. It has an envelope delay characteristic given by 

D 3 (co) = (D m /2)[1 + cos a>772] for | « | < 2t/T (126) 

and hence a phase characteristic 

3 ( w ) = (D m /T)[(a>T/2) + sin uT/2] for | co | < 2w-/T. (127) 

The three envelope delay chacteristics are shown in Fig. 19. 

The general expressions for the error rate, derived previously for the 
various modulation systems, have been evaluated on a digital computer, 
with the equivalent receiving filter having each of the above phase char- 
acteristics. Values of maximum delay up to 10 symbol durations have 
been considered. The resulting error rates are compared with those 
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Fig. 19 — Envelope delay characteristics. 
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when no delay is included, and the improvement is plotted, as a function 
of the maximum delay, in Figs. 20 through 23. 

For all but the differentially coherent phase shift keying systems, the 
delay networks have no effect for a value of a equal to 1. This has been 
discussed in the preceding sections; by Parseval's theorem the phase 
characteristics of the noise have no effect when a is equal to 1. As the 
value of a increases, the networks become more effective. This phenome- 
non is to be expected. When the noise burst is spread over many bits, 
the peak value is reduced and many bursts which would have caused 
errors no longer do so. However, this improvement is reduced somewhat 
by the fact that bursts of very large amplitude remain large enough, 
even after spreading, to cause errors. Such bursts, which would have 
caused one or two errors, now cause many. As a increases, the percentage 
of such high amplitude bursts decreases, and the improvement increases. 
This effect, which reduces the over-all improvement to be gained by the 
use of complementary delay niters, prompts the consideration of limit- 
ing at the input to the delay filter at the receiver. If the limiter is set to a 
value just above the peak value of the signal, bursts of very high ampli- 
tude are clipped and after spreading, will not cause as many errors as 
they would have without limiting. 

A precise analysis of the effect of limiting prior to spreading is not 
practicable. The process is nonlinear, and the response of the combined 
limiter-filter depends strongly on the amplitude K and time of occurrence 
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Fig. 20 — Effect of complementary delay networks on double sideband AM 
systems. 
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Fig. 21 — Effect of complementary delay networks on single-sideband AM 
system. 

r of the noise burst, as well as on the particular data sequence being 
transmitted. In addition, the results of such an analysis would be quite 
sensitive to the model chosen to represent the channel. In the systems 
discussed in this paper, the transmission medium and receiving filter 
have both been linear, and it was therefore possible to combine their 
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Fig. 22 — Effect of complementary delay networks on frequency shift keyed 
system. 
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Fig. 23 — Effect of complementary delay networks on differentially coherent 
phase shift keyed system. 

characteristics for the purpose of analysis. The channel was considered 
to have unity gain in the band of frequencies passed by the receiver 
band-pass filter; the transmission characteristics of the channel at fre- 
quencies outside that band had no effect on the signal or noise at the 
output of the filter. If limiting is introduced, however, it is most effective 
if it is done at the point where the noise bursts are most impulsive. This 
point is at the input to the receiver, prior to any filtering or spreading, 
where their spectrum is widest. With a nonlinear device between the 
channel and the receiving filter, their characteristics cannot be com- 
bined. The effect of the limiting depends very much on the transmission 
characteristics of the channel outside the pass band of the filter. If the 
channel has a considerably wider band than the filter, the limiting will be 
much more effective than if the bandwidths are comparable. 

For these reasons, an analytic evaluation of the improvement result- 
ing from the introduction of limiting is not practicable. An evaluation 
by simulation techniques or experimental procedures is more feasible. 
Limiting is discussed in this paper only to point out that further improve- 
ment is possible. 
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